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ABSTRACT 

The disk around AU Microscopii, an Ml star in the (3 Pictoris Moving Group, is extraordinarily well- 
suited for comparison to the (3 Pic debris disk (type A5V). We use far-UV absorption spectroscopy of 
AU Mic to probe its edge-on disk for small amounts of H2, the primary constituent of gas giant planets. 
Our conservative upper limit on the line-of-sight H2 column density is 1.7 x 10 19 cm" 2 , which is 18.5 
times lower than the limit obtained from non-detection of sub- mm CO emission (Liu et al. 2004). In 
addition, there is a hint of H2 absorption at a column density an order of magnitude or more below 
our upper limit. The RVto-dust ratio in the AU Mic disk is < 6:1, similar to that in the (3 Pic disk. 
This shows that the primordial gas has largely been dissipated in less than about 12 Myr for both disks, 
despite their very different stellar masses and luminosities. It is extremely difficult to form a giant planet 
around AU Mic with current core-accretion scenarios in such a short time. 

Subject headings: planetary systems: protoplanetary disks, stars: low-mass, stars: individual (AU 
Microscopii) 
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disk (Kalas et al. 2004). This disk is uniquely suited for 
comparison to the (3 Pic debris disk and provides a valu- 
able opportunity to examine the effect of stellar mass on 
disk evolution. 

The inclination of the AU Mic disk within 50 AU of the 
star is less than 1° (Krist et al. 2005), even more edge-on 
than the (3 Pic disk (i = 3°; Heap et al. 2000). CS gas in 
the (3 Pic disk has been observed with line-of-sight absorp- 
tion spectroscopy (e.g. CO; Roberge et al. 2000), and a low 
upper limit on the H2 abundance was set using far-UV ab- 
sorption spectroscopy (Lecavelier des Etangs et al. 2001). 
This limit showed that the primordial gas has largely been 
dissipated. In this paper, we apply similar analysis to far- 
UV spectra of AU Mic, to probe sensitively for traces of 
primordial H2 gas. 

2. OBSERVATIONS 

We used two sets of archival far-UV spectra of AU Mic, 
one each from the Far Ultraviolet Spectroscopic Explorer 
(FUSE) and the Hubble Space Telescope Space Telescope 
Imaging Spectrograph (STIS). Detailed descriptions of the 
data and analyses of the stellar emission lines may be 
found in Redfield et al. (2002) and Pagano et al. (2000). 
AU Mic was observed with FUSE on 2000 August 26 
and 2001 October 10 for 17.3 ksec and 26.5 ksec, respec- 
tively. These spectra cover 905 - 1187 A with a resolution 
of about 15 km s" 1 (R = A/AA = 20,000 ± 2000). The 
data were recalibrated with CALFUSE version 2.4.0. The 
STIS data were taken 1998 September 6 (exposure time 
= 10.1 ksec). These spectra cover 1170 - 1730 A with a 
resolution of 4.4 km s" 1 at 1350 A (R = 68,000). We 
recalibrated the data with CALSTIS 2.17b. 

Stellar flares occurred in both datasets. Since emission 
line shapes and strengths can change during flares (e.g. 
Redfield et al. 2002), data taken at these times were ex- 
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1. INTRODUCTION 

A circumstcllar (CS) disk around a young star evolves 
from a massive remnant of star formation (a primordial 
disk, composed mostly of H2 gas) to a relatively low mass, 
gas-free planetary system. But there are difficulties rec- 
onciling observations of CS disks with current theories of 
giant planet formation. In the traditional picture, giant 
planets form by accretion of gas onto a massive solid core, 
which must form before the primordial gas in the disk 
dissipates. Core formation takes at least a few Myr for 
solar-type stars, much longer for lower mass stars (Laugh- 
lin et al. 2004). Observations suggest that the primordial 
gas dissipation timescale is shorter than about 10 Myr, 
possibly as short as 1 Myr in regions of high- mass star for- 
mation (Bally et al. 1998). With these short gas lifetimes, 
it would be fairly difficult to form giant planets with cur- 
rent core-accretion scenarios around most solar-type stars 
and extremely difficult around all low-mass stars. 

However, giant planets have been found around many 
G and K-type stars, and two M-type stars (e.g. Butler 
et al. 2004) . This has led some workers to suggest alterna- 
tive theories for giant planet formation, i.e. that they form 
rapidly by direct collapse in gravitationally unstable disks 
(Boss 1998). But it is not clear that current core-accretion 
models must be rejected, due to the uncertain timescale for 
primordial gas dissipation around stars of various masses 
in different star- forming environments. Comparison be- 
tween the abundance of giant planets around M-stars and 
the typical H2 gas lifetime in their disks will be a powerful 
constraint on theories of giant planet formation. 

The Ml star AU Microscopii (GJ 803) has long been 
studied as a nearby flare star (d = 9.94 ± 0.13 pc). It is 
co-moving with the A5V star (3 Pictoris, indicating they 
formed in the same region and are the same age (12^4 Myr; 
Zuckerman et al. 2001). The CS dust around AU Mic was 
imaged in reflected light and shown to lie in an edge-on 
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eluded from our analysis. The FUSE bandpass is cov- 
ered by 8 partially overlapping spectra, two of which have 
the same zero-point wavelength offset (LiFIA and LiFIB). 
The bandpasses of the LiFIB FUSE spectrum and the 
STIS spectrum overlap, so we set the absolute wavelength 
calibration of the LiFIA and LiFIB spectra by cross- 
correlating on the C III A1176 stellar emission line seen 
in both datasets. The zero-point offset of the FUSE spec- 
tra was small (w 5 km s _1 ). 

3. ANALYSIS 

Some of the many strong far-UV electronic transitions 
of H2 overlap with commonly seen stellar emission lines, 
i.e. the O VI AA1032, 1038 doublet originating in the stellar 
transition zone, and the chromospheric C II AA1036, 1037 
doublet. The portion of the FUSE LiFIA spectrum includ- 
ing these lines is shown in Figure 1 . We have demonstrated 
that these emission lines can provide the background flux 
for absorption spectroscopy of H 2 (Roberge et al. 2001; 
Lecavelier des Etangs et al. 2001). 

3.1. Creation of the Emission Line Model 

The first step in analyzing the FUSE spectra for molec- 
ular hydrogen absorption was to model the intrinsic flux in 
the emission lines. This is an emission model with infinite 
spectral resolution that will be multiplied by an intrinsic 
absorption model, then convolved with the instrumental 
line-spread function (LSF) to produce a model of the ob- 
served flux. 

The O VI A1032 line is little affected by H 2 absorp- 
tion, and the doublet lines are always in the optically thin 
ratio of 2:1, unless the A 1038 line is suppressed by H2. 
The measured integrated flux ratio of the two O VI lines 
is 1.97 ± 0.19, immediately indicating that there is little 
line-of-sight H 2 absorption. We began by least-squares fit- 
ting the sum of a narrow and a broad Gaussian to the 
A1032 line, since chromospheric and transition zone emis- 
sion lines from an active star like AU Mic are not well 
fit with single Gaussians (e.g. Redfield et al. 2002). The 
widths and peak intensities of the Gaussians were adjusted 
to remove the effect of the Gaussian FUSE LSF. The pa- 
rameters from the fit to the A 1032 line were used to model 
the A1038 line, with the peak flux set by the optically thin 
ratio. 

The C II AA1036, 1037 lines provide a tighter constraint 
on H2 absorption, but were difficult to model, since both 
can be affected. Therefore, we analyzed the C II A1335 
triplet in the STIS data, shown in Figure 2. Examination 
of the line profiles shows that the blue wing of the A1334.5 
emission line is cut off by ground-state C II absorption, 
which was initially identified as interstellar (IS), but the 
A1335.7 emission line is not affected by excited-state C II 
absorption. 

We found the parameters of the C II emission by analysis 
of the A1335.7 line, which is blended with the weak A1335.6 
line. The emission was modeled as the sum of a narrow and 
broad Gaussian for each line, setting the line separation 
from the rest wavelengths. The relative line strengths were 
set using the products of the upper state multiplicities and 
Einstein A-coefficients. We then performed \ 2 minimiza- 
tion. The widths and peak intensities from this analysis 
were adjusted to remove the effect of the STIS LSF, which 



is very small since the lines are well-resolved. The narrow 
component is at v n = —3.0 ± 3.4 km s -1 , where the error 
bar is the ±la accuracy of the STIS wavelength calibration 
after a standard target acquisition (Pagano et al. 2000). 
Within the uncertainty, the C II emission is at the velocity 
of the star (v* = — 4.89±0.02 km s _1 ; Barrado y Navascues 
et al. 1999). As is usually the case, the broad component 
is slightly redshifted (v^ = 3.0 ±3.4 km s _1 ). The narrow 
component has FWHM = 35.63 ± 0.55 km s" 1 and the 
broad component has FWHM = 129.4 ± 3.9 km s _1 . 

The parameters from the A1335.7 minimization were 
used to model the A1334.5 emission line, again setting the 
wavelength separation and relative line strength from the 
atomic data. We found that the A1335.7 line is not opti- 
cally thin, since the peak intensities of the Gaussians in 
the A1334.5 line had to be increased slightly to fit the blue 
wing of the line. The emission model of the triplet con- 
volved with the STIS LSF is shown with a dashed red line 
in Figure 2. 

To find the parameters of the C II A1334.5 absorp- 
tion, we used a Voigt intrinsic line profile to create an 
absorption model. A total model was made by multiply- 
ing the triplet emission model by the absorption model, 
then convolving the result with the STIS LSF. The ab- 
sorption parameters from x 2 minimization of the total 
model are N = x 10 14 cm~ 2 (column density), 

b = 17.50 ± 0.63 km s _1 (Doppler broadening parameter), 
and v = — 39.5±3.4 km s _1 . The convolved final model for 
the C II A 1335 triplet is overplotted with a solid red line 
in Figure 2. The C II absorption is broad and blueshiftcd 
with respect to both the star and other absorption lines 
seen toward AU Mic (e.g. D I at v = —21.7 km s _1 ; Red- 
field & Linsky 2004). This suggests that the line-of-sight 
C II absorption may not be primarily IS. 

Next, we modeled the C II emission doublet in the FUSE 
spectrum, using the velocities and widths of the narrow 
and broad Gaussians and the ratio of their peak intensi- 
ties from the STIS analysis. The separation of the two 
lines and their relative strengths were calculated from the 
atomic data. The doublet emission model was multiplied 
by a C II A1036 absorption line, using the parameters given 
in the previous paragraph. The remaining free parameter 
was a scaling factor on the peak intensity of the A1036 
line, which was determined by matching the model to the 
blue wing of the line. We found that the ratio of the line 
strengths calculated from the atomic data had to be ad- 
justed slightly to fit the blue wing of the A1037 line, again 
indicating it is not optically thin. 

The C II and O VI models were added together to pro- 
duce an infinite resolution emission model for the 1030 - 
1040 A region of the FUSE spectrum. This model con- 
volved with the FUSE LSF is plotted in Figure 1 with 
a solid blue line. One can immediately see that the C II 
emission lines in the FUSE data arc slightly narrower than 
expected from the width of the STIS C II lines. There is 
no intrinsic reason for this, since we excluded flares from 
both data sets. Self- absorption suppresses the peak flux 
of an emission line, but does not make it appear narrower. 
The emission lines are resolved, so the uncertainty in the 
spectral resolution of the FUSE data has little effect on 
the model. Other obvious FUSE instrumental or data ar- 
tifacts tend to make emission lines appear broader than 
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they really are, not narrower. It seems there are only two 
ways to make the FUSE C II lines appear narrower than 
the STIS lines. The first is a calibration error that artifi- 
cially broadened the STIS lines. There is such an error in 
STIS far-UV time-tag data (like ours) retrieved from the 
archive before 2004 April 21 4 . However, we recalibrated 
our data with the corrected pipeline software. The other 
way is to remove some emission from the FUSE lines with 
superimposed H2 absorption (discussed further in §4.2). 

3.2. Molecular Hydrogen Model Fitting 

An H2 absorption model was created using Voigt line 
profiles and molecular data from Abgrall et al. (1993). 
The free parameters were N(J = 0), the column den- 
sity in the v — 0, J = level, N(J = 1), the column 
density in the v = 0, J = 1 level, and v, the velocity of 
the gas. Only the two lowest energy levels were included, 
since they contain more than 90% of the molecules at the 
typical temperatures of debris disks (< 200 K). The in- 
trinsic width of any absorption is unknown, so we made 
the most conservative assumption (unresolved lines) and 
set the Doppler broadening parameter to be very small 
(6 = 1 km s _1 ). The H 2 model was multiplied by the 
FUSE emission model, then convolved with the FUSE 
LSF. The resulting total model was compared to the data 
between 1035.5 A and 1037.7 A (the region most affected 
by H2 absorption). To ensure that all local minima were 
found, x 2 minimization was performed using very large 
ranges of parameters (10 12 cm~ 2 < N(J = 0) and N(J = 
1) < 10 22 cm~ 2 , -60 km s" 1 < v < +60 km s _1 ). 

4. RESULTS 
4.1. Upper Limits on H2 

Using our contours of x 2 , we set upper limits on the 
line-of-sight H 2 column densities of N(J = 0) < 1.6 x 
10 19 cm~ 2 and N(J = 1) < 6.3 x 10 17 cm" 2 . These col- 
umn densities are ruled out at the > 6a level for all v. In 
Figure 3, the convolved H2 absorption model with these 
values is shown, and the total model overplotted on the 
data near the C II lines. At lower column densities, the 
contours of x 2 are complicated by the presence of a lo- 
cal minimum in addition to the absolute minimum, which 
is why we set the limits more conservatively from the 6a 
contour rather than the typical 3a contour. The upper 
limit on the total line-of-sight H2 column density toward 
AU Mic is Atotai < 1-7 x 10 19 cm~ 2 ; we did not assume 
a gas temperature to get this limit. If one assumes that 
the H2 energy levels are thermally populated and fixes the 
gas temperature at any value below a few thousand K, the 
total upper limit decreases. 

From sub-mm observations, Liu et al. (2004) placed a 
3a upper limit on the CO column density in the AU Mic 
disk of 6.3 x 10 13 cm" 2 , assuming any gas was in thermal 
equilibrium with the observed dust (T = 40 K). Using 
a CO/H2 ratio of 10~ 7 , this corresponded to a total H 2 
column density of 6.3 x 10 20 cm" 2 . The CO/H 2 value 
is unusually low compared to the typical IS value, since 
the authors tried to take into account the fact that the 
AU Mic disk is optically thin to IS UV photons. To con- 
vert this column density from emission observations to a 

4 see http://www.stsci.edu/hst/stis/documents/newsletters/stan0412. 



line-of-sight absorption column, we divided the sub-mm 
upper limit by 2. This gives the lowest possible value, by 
assuming that the radial density distribution of the gas is 
extremely centrally peaked. Our x 2 contours show that 
the sub-mm upper limit on the line-of-sight H2 column 
density is ruled out at the > 20ct level for all v. Such 
a large column density would dramatically suppress the 
O VI A1038 emission line and completely obliterate both 
C II emission lines (see Figure 1). Our upper limit in the 
previous paragraph is 18.5 times lower than the sub-mm 
upper limit, docs not assume a gas temperature, and is 
free from the large uncertainties introduced by assuming 
a CO/H2 ratio. 

4.2. Possible Line-of-Sight H2 Absorption 

The fact that the FUSE C II emission lines are narrower 
than the STIS lines suggests there is a small amount of H 2 
absorption toward AU Mic. The minimization with fixed b 
described above showed a second minimum in the x 2 con- 
tours at H2 column densities about an order of magnitude 
below our 6a upper limits. In addition, there have been 
claims of weak fluorescent H2 emission lines in the FUSE 
and STIS spectra (Pagano et al. 2000; Redfield et al. 2002); 
however, these lines are only detected at the « 2a level. 

We attempted to detect line-of-sight H 2 absorption by 
performing x 2 minimization with all parameters free, in- 
cluding b. These contours showed a minimum at large 
b and column densities a few orders of magnitude below 
the upper limits in §4.1. But as before, there were multi- 
ple minima in the x 2 contours, indicating there are multi- 
ple sets of parameters producing H2 models which fit the 
data about equally well. Therefore, the model is overde- 
tcrmincd given the quality of the data. But taken at face 
value, our analysis indicates that the putative line-of-sight 
H 2 absorption is broad and blueshifted with respect to the 
star, as was the line-of-sight C II absorption in the STIS 
data. 

5. DISCUSSION 

Liu et al. (2004) measured a sub-mm dust mass of 
0.011 M e within 70 AU of AU Mic. Their upper limit 
on the H2 mass (1.3 M®) gives a limit on the H2-to-dust 
ratio of < 118:1. Our limit on the line-of-sight H2 column 
density reduces the possible gas mass by at least a fac- 
tor of 18.5, giving an upper limit on the H 2 -to-dust ratio 
of < 6:1. The typical IS ratio is 100:1; therefore, the H2 
gas in the AU Mic disk has been depleted relative to the 
dust. The primordial gas is largely dissipated, and the gas 
lifetime was less than about 12 Myr. 

Laughlin et al. (2004) showed that a core massive 
enough to gravitationally accrete a large gaseous envelope 
could not form within 10 Myr around a 0.4 M Q star. The 
only way current core-accretion scenarios can form a gi- 
ant planet around an M-star is if its primordial gas sur- 
vives much longer than about 10 Myr. Therefore, it is 
unlikely that a giant planet could form around AU Mic 
through these scenarios. There is indirect evidence of a gi- 
ant planet around AU Mic, since the central region of the 
disk is cleared of dust grains (Liu et al. 2004) and the disk 
appears slightly warped in the HST-ACS scattered light 
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image (Krist et al. 2005). Central holes and warps have 
been attributed to the dynamical influence of an unseen 
giant planet. 

We now consider the case of (3 Pic. The disk dust mass 
is 0.04 M® (Dent et al. 2000). Using the upper limit on 
the H 2 gas mass from FUSE (< 0.1 M e ; Lecavelier des 
Etangs et al. 2001), the (3 Pic H2-to-dust ratio is < 3:1, 
similar to the limit for AU Mic. This indicates that the 
primordial gas lifetimes in the two systems were both quite 
short, despite the very different stellar masses and lumi- 
nosities (I/au Mic = 0.1 L Q ; L/jpic = 8.7 L Q ). There is 
no O or B star in the (3 Pic Moving Group, so the short 
lifetimes are not easily explained by external irradiation. 
Given their spectral types, photoevaporation by the cen- 



tral stars seems unlikely to produce short gas lifetimes in 
both systems. Other possible dissipation mechanisms in- 
clude stellar winds, magnetospheric accretion, and planet 
formation. 

Our observation of broad, blueshifted line-of-sight C II 
absorption (and possibly H 2 ) would lead one to suspect 
that there is outflowing, low-ionization gas near the star, 
possibly in a stellar or disk wind. However, there may be 
an unusually large number of blended IS components at 
different redshifts along the sight line to AU Mic. We plan 
to investigate further by re-examining the putative fluores- 
cent H2 emission lines and all atomic absorption features 
in the STIS spectrum previously identified as IS. 
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Fig. 1. — Portion of the FUSE LiFIA spectrum of AU Mic with a demonstration of the effect of line-of-sight H2 absorption on the O VI 
AA1032, 1038 and C II AA1036, 1037 emission doublets. The intrinsic emission model, convolved with the FUSE LSF, is ovcrplotted with 
a blue solid line. A convolved H2 absorption model with the parameters of the upper limit from non-detection of sub-mm CO emission 
(7V tota j = 3.15 X 10 20 cm -2 , T = 40 K; Liu et al. 2004) is plotted with a black dashed line. For the H2 model shown, we assumed unresolved 
lines and set the velocity of the gas to the stellar velocity (v+ = —4.89 km s _1 ). The total model is the absorption model multiplied by the 
emission model, then convolved with the FUSE LSF (red solid line). The O VI A1032 emission line is little affected by H2 absorption, but 
such a large H2 column density would dramatically suppress the O VI A1038 line and completely absorb both C II emission lines. 
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Fig. 2. — Portion of the E140M STIS spectrum of AU Mic showing the C II A1335 emission triplet (analyzed to aid in modeling of the 
FUSE C II emission lines). The blue wing of the C II A1334.5 line is cut off by linc-of-sight C II absorption. The narrow and broad Gaussians 
from our modeling of the C II A1335.7 line are overplotted with dotted blue lines. The emission line model of the triplet, convolved with the 
STIS LSF, is overplotted with a red dashed line. The convolved C II A1334.5 absorption model is plotted with a black dashed line. The total 
model is the absorption model multiplied by the emission model, then convolved with the STIS LSF (solid red line). 
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Fig. 3. — Upper limit on line-of-sight H2 absorption toward AU Mic. The portion of the FUSE LiFIA spectrum near the C II AA1036, 1037 
doublet is shown, with the convolved emission model overplotted (blue solid line). The black dashed line shows the convolved upper limit H2 
absorption model (N(J = 0) = 1.6 X 10 19 cm -2 , N(J = 1) = 6.3 X 10 17 cm -2 ). The total model is the absorption model multiplied by the 
emission model, then convolved with the FUSE LSF (solid red line). For the H2 model shown, the velocity was set to +4.0 km s . This is 
the most redshifted velocity expected for line-of-sight CS gas in Keplerian rotation in the edge-on disk, given the radial velocity of the star 
(v* = —4.98 km s _1 ) and the combined accuracy of the FUSE and STIS wavelength calibrations (±8.7 km s _1 ). H2 absorption models with 
velocities blueward of +4.0 km s — 1 have an even more dramatic effect on the C II emission lines. 



